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Abstract 



We propose a simple friction model for isolated polymer chains on a solid 
substrate. The chains are pulled at constant velocity by one end, the other 
end can be trapped on the solid substrate on localised sites. We focus on the 
energy dissipation due to the traps. This simple model leads to non trivial 
friction laws, depending on the velocity and the distance between traps. Some 
refinements of the model such as the effect of thermal fluctuations are also 
reported. 
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Friction phenomena occur at the interface between two media under relative motion. 
Their understanding is crucial for many technological applications and raises many funda- 
mental questions. This is one of the oldest problems in physics but there is still a relatively 
poor understanding of the microscopic mechanisms of friction. 

Dry solid friction is the older known form of friction, first introduced by da Vinci and 
Amonton. The dry friction laws give the friction force in terms of the load and the static and 
dynamic friction coefficients /i s and fid- A first surprising result is that the friction coefficients 
do not depend of the contact area. The first microscopic model of Amonton's law |l],[2j 
appeared in the middle of this century and shows the importance of the interface roughness. 
The dynamic friction coefficient depends only very weakly on the sliding velocity ti |-|. 
and in practice, the friction force reaches a constant value at vanishing velocity. There is 
still a poor understanding of dry friction, both from the theoretical and the experimental 
points of view, in part due to the difficulties to control the interfacial properties @. 

Fluid friction can be studied within the framework of hydrodynamics and is better un- 
derstood. For a solid moving in a newtonian fluid, the friction law involves the viscosity 
of the fluid rj, its density p, its size R and the velocity u. The dimensionless Reynolds 
number R e = gives the ratio between the inertial and viscous forces. A simple scaling 
argument leads to the force acting on the solid as: / ~ r]Rug(R e ), where g is a function of 
the Reynolds number. At low velocities, ie low Reynolds numbers, the friction force must 
be independent of the density p, and thus is proportional to rjRu, (Stokes friction). At high 
Reynolds number, the viscous stresses are negligible and / ~ rjRuR e ~ pR 2 u 2 . 

The friction between two polymeric media is a subject of major interest for many prac- 
tical applications (lubrification, adhesion, extrusion, wetting). Polymeric systems also allow 
to work on well controlled surfaces and several microscopic parameters that can influence 
friction can be varied, such as a chain grafting density, chain lengths, chain affinity ... On 
the other hand, the complexity of the dynamic of the interfacial chains leads to a large 
variety of frictional behaviors: gel or me lt friction, slip transition |T0|-|I2]1 and stick-slip 

MM- 

In a recent paper |TjJ, Subbotin et al proposed a model describing the friction of confined 
polymers melt. In their model, the energy dissipation is simply due to the high friction of 
monomer anchored on the surface and the effect of loops and bridges on the friction law were 
precisely studied. In this paper, we present a simple friction model for a single polymer chain 
on a solid surface, focusing mainly on the microscopic mechanism of energy dissipation. We 
consider a polymer chain in a solvent, pulled by one of its extremities (s=N where s is the 
monomer number) at a constant velocity u (figure |l|). The other extremity (s=0) can be 
trapped on the surface by localised defects. Our aim is to calculate the dissipated energy 
due to the traps. This extremely simple and naive model gives a non trivial behavior of the 
friction force which is similar to solid friction in the sense that it may not vanish in the limit 
of small velocities. It could be also possible to test this behavior by micromanipulation of 
isolated polymer chains in a solvent with an atomic force microscope. Finally, it could be a 
possible way to introduce an energy dissipation in molecular motor models [16[ . 

We first present the one chain problem on a linear and periodic array of traps. In this 
first part, we simply assume that the trapped monomer is leaving the trap when the tension 
reaches a critical value r. In the second part, we show how thermal fluctuations can induce 
the desorption of this monomer trapped in a potential barrier. 
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In a first approximation, the dynamics of a polymer chain can be described by the well- 
known Rouse model [|l7j], which neglects hydrodynamic interactions. The polymer chain is 
represented as a set of N beads at position z connected by springs with entropic rigidity 
k = ^r~, where a is the monomer size. Each monomer behaves as a solid particle with a 
friction coefficient ( proportional to the solvent viscosity 77, £ ~ rja. The Rouse equation 
which describes the dynamics of the chain is written(in one dimension) as: 

C*M-*£.M + /M (i) 

where z(s,t) is the position of monomer s at time t and f(t) is a Langevin random 
force. The slowest relaxation mode of the chain has a characteristic time (the Rouse time) 
Tr = . As it ignores hydrodynamic interactions, the Rouse model does not give a good 
description of polymers in dilute solutions. Nevertheless, the Rouse description is concep- 
tually important and we choose it as a starting point to describe our results on polymer 
friction. In appendix A we show how it is possible to take into account hydrodynamic in- 
teractions within the framework of the Zimm model. This complication does not change 
qualitatively the behavior of the friction force. Our results are obtained for gaussian chains 
in a G solvent. They can easily be generalised to polymer in a good solvent, but for a sake 
of clarity we do not show the results in the present paper. 

Our paper is organised as follows. In the first section, we consider the interaction between 
one chain and a linear array of traps along the direction of the velocity, ignoring any thermal 
fluctuations. The traps are pointlike on the surface and the chain end point remains trapped 
as long as the chain tension is smaller than a critical value r which is the maximum force 
that can be exerted on the chain end by the trap. We first discuss the single trap problem 
and then a periodic one-dimensional array of traps with a period d. In appendix D we relate 
this model to the more realistic problem where the surface has a finite density of traps n per 
unit area. Complete calculations can be performed within the Rouse model; in the main text 
of the paper, we rather give simple scaling arguments postponing the detailed calculations 
in an appendix. Section [TT] is devoted to the effect of thermal fluctuations. We describe 
each trap more realistically as a potential well and we estimate the escape probability of 
the trapped chain end when the tension increases. This amounts to calculating an effective 
critical tension with which the results of the previous section can be used. Following reference 
|18| , we describe the release of the chain end as a Kramer process and we use a self-consistent 
force approximation. The last section summarises our results and discusses some possible 
extensions. 

I. POLYMER CHAIN INTERACTING WITH TRAPS 

A. The one trap problem. 

We study a Gaussian polymer chain interacting with one trap on a solid surface. For 
simplicity, we consider the chain as two-dimensional and lying in the surface with a Gaussian 
conformation (at rest). The end point of the chain s = N moves at a velocity u. When 
the chain is free (not trapped), its average conformation can be obtained from the Rouse 
equation |I[ The average position Zf ree (s,t) of monomer s at time t is 
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l(u 2 l fs\ 2 

Zfree{S,t) = -—3 = ~uTr l^ — J (2) 

The tension of the chain at monomer s is given by the friction force on the chain end 
section containing s monomers r c h(s, t) = k dz f T ^ s > t > = (us. The tension on the free end of 
the chain vanishes, and the external force applied to pull the chain is N(u = j^uTr. This 
can be viewed either as the friction force on the N monomers, or as the elastic tension of 
N springs in series extended over a distance uTr. In the absence of trapping, the energy 
dissipated per unit time by the solvent viscosity is simply = N(u 2 . 

When the end monomer (s=0) gets trapped, the chain is progressively elongated by the 
motion of the pulled end (s=N). The tension applied on the trapped monomer increases 
until it reaches the critical value r. The trapped chain end is then released and the chain 
contracts back to its free conformation. 

The Rouse equation for the trapped chain allows the calculation of the chain tension 
(appendix : 



Tch{s,t) = k^- 



t s 2+^1/ 2 2 t A / n 



(3) 



The tensions on the trapped end (r tr (t) = r c h (s = 0, t)) and on the pulled end 
{j c h (s = N,t)) are plotted on figure 0. At short time after the trapping (t <^ Tr), the 
tension on the pulled end remains constant (equal to the viscous friction force N(u), and 

/ \ 1/2 

the tension on the trapped monomer increases as r tr ~ (Nu ( ^ ) • At larger times, 
(t ^> Tr), the tension is uniform along the chain and both tensions increase as (Nu^. 
After the sticking time t s , the tension r tr becomes equal to the critical value r. 

The characteristic velocity of the problem is obtained by comparing the friction force on 
the chain (Nu to the maximum force r that can be exerted by the trap u c = = ^p^. 

If u u c , the viscous force on the chain is smaller than r. When the trapped end is 
released from the trap, the whole chain is stretched with a constant tension and its extension 
is L = The tension acting on the first monomer at time t is r tr (t) ~ j?ut and the sticking 
time t s is given by : 

t s = —Tr (U < U c <^> t s > Tr) 

u 

The work performed to extend the chain is After the release of the trapped end, 
the elastic energy stored in the chain Wd = ^f- is dissipated by a viscous process. This 
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argument is essentially identical to the classical argument of Lake and Thomas |19|] which 
has been used to study the adhesive properties of polymers by de Gennes [p0| , pTt1 . 



If u 3> u c , the viscous force on the chain is larger than r, and only a small part of the 
chain involving So(t) ~ monomers is elongated when the release of the trapped end 

occurs. It contributes to the tension on the trapped monomer as r tr {t) ~ j^fj u ^ anc ^ ^ ne 
sticking time t s is now given by : 
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where the numerical prefactor is obtained from the detailed calculations of the appendix. The 
force acting on the pulled monomer remains constant and equal to N(u during the trapping 
with a good approximation. After the release of the traped end, the tension of the pulled 
monomer is calculated in appendix |B] and is given by ( p4|) . The extra energy dissipated 
due to the trapping is of the order of the elastic energy stored during the trapping in the 
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deformed part of the chain Wd = So(i s )-T- = N 1 ^—. More complete calculations confirming 
these scaling law can be done using the Rouse model. They are presented in appendices |A] 
and |B|. After the release of the anchored monomer, the chain relaxes to the steady state, 
and its elastic energy is dissipated. The total dissipated energy is 



We now consider a polymer pulled at a velocity u on a two dimensional surface with a 
density of traps n. We refer to appendix [D] for a discussion of the trapping probability that 
gives the mean distance d(u) between trapping events. The problem is then equivalent to a 
one dimensional array of traps along the velocity direction with an average density l/d(u). 
For simplicity, we consider here a periodic lattice with a distance d(u) between traps. We 
then calculate the average friction force in a steady state as Ff r = ^ where Wd is the 
energy dissipated on each trap. 

We first consider the limit of very low velocities u <C u c . When the chain is trapped 
on one defect and is just on the verge of detaching, the chain has a constant tension r 
and its size is L = If this is smaller than the distance d between traps, the traps are 
independent and we can use the results of the previous section. The energy dissipated per 
trap is Wd ~ ^r- and the friction force is Ff r ~ If d <C Nr/k, the chain end gets 

trapped before the full relaxation of the tension. The number of monomers where the tension 
has relaxed increases as t 1 / 2 and the position of the first monomer is given by z{t) ~ 7^172 • 

The chain gets trapped again when z(t) = d or after a time t a ~ • The tension then 
becomes uniform along the chain over the Rouse time Tr and equal to r' = k(L — d)/N. 
The chain stretches then back to a uniform tension r over the time d/u. The dissipated 
energy per trap is Wd = ^ and the friction force is Ff r ~ r. 

This quasi-static description remains valid as long as d/u > Tr. If this is not the case, 
the tension does not have time to relax to the smaller value r' over the whole chain and 
remains equal to r (in fact sightly larger) on the pulled monomer. As long as the viscous 
force is small (u < u c ) this tension balances the friction force and thus Ff r ~ r. In this case, 
the dissipated energy per trap is of the order of the elastic energy of a chain spanning the 
distance d i.e. containing s = kd/r monomers. 

If the velocity u is larger than the critical value u c , the total dissipation is dominated 
by the viscous forces. The defects are independent if t 8 (u) <C ^ i.e. if d > In this case 

the dissipation is due to individual traps and the friction force is Ff r ~ ~ 1ST if- ^ 

2 

decreases with the velocity. If d tV, as above, the dissipated energy per trap is of the 




(4) 



B. The many-traps problem. 
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order of the elastic energy of the chain section spanning the length d with a tension r and 
the friction force is Ff r = t. 

All these results are summarized on figure |3] in a velocity-distance between traps (u, d) 
diagram. 



II. THERMAL FLUCTUATIONS. 

In the previous section, we have considered point-like traps and we have ignored thermal 
fluctuations describing the chain statistics by the average position of each monomer. In a 
more realistic model, these two assumptions are relaxed. Each trap is an attractive potential 
well V(z — zq) for the chain end with an extension b around the trap position z and a depth 
Uq. The maximum tension that can be exerted by the trap r is the slope of this potential 
at the inflection point Z\. We want to study here strong trapping and we only focus on the 
case where 1. This can be achieved for a reasonable tension if the size of the trap is 

larger than the monomer size a (typically a few times) but is much smaller than the chain 
size and the distance between traps d. A typical trapping potential is sketched on figure f|.a. 

Thermal fluctuations have two effects; they can induce a new trapping of a chain end 
that has just escaped from the trap; they can also cause the release of the end monomer 
even if the tension is lower than the critical value r. 



In references [p2jP3| , Wittmer et al and Clement et al have shown that the readsorption 
of the chain can occur until the terminal monomer leaves the first Pincus blob of size £ = 
ksT /r tr - The size of this Pincus blob in the strong trapping limit is smaller than the size of 
the trap b when the tension is equal to the critical tension. The retrapping effect can then 
be neglected. 

The escape of a particle submitted to an external force from a potential well can be 
studied using Kramers rate theory. If the constant external force is r tr , particle is submitted 
to the effective potential U(z) = V(z) — r tr z (choosing the origin at the center of the trap). 
This potential has a minimum value at a position z min corresponding to the equilibrium 
position of the particle and a maximum value at a position z max (z min < Z\ < z max ). The 
escape time due to thermal fluctuations from the potential is the first passage time of the 



particle at the position z max . Within Kramers rate theory pl| , this is given by 

t e {rtr) = tiexp \{U(z max ) - U(z min )} /k B T] (5) 
where t\ is given by the curvature of the potential U(z) at the positions z m i n and z max 

*iW = " N i/ 2 (6) 

d 2 U d 2 U \ 

TP TP ; 

Our aim here is to use these result to discuss the escape of the trapped end under the 
action of the time-dependent tension and to calculate the chain tension r e when the trapped 
end escapes from the potential well. The results of the previous section can then be used 
with the effective critical tension r e instead of r. 

In reference two approximations (so-called adiabatic and self-consistent force ap- 
proximation) are proposed to discuss the escape of a particle from a potential well under 
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the action of a non constant force. We have checked that both lead to very similar results. 
For a sake of clarity, we present here the simpler self-consistent force approximation. In this 
approximation, the tension is obtained self-consistently by calculating the escape time from 
equation [| and |6| for a tension r tr increasing with time and by imposing that r e = r ir (t e ). 

In the limit of small velocities u <C u c , the chain tension at the end point increases as 
Ttr = nrr- The self-consistency equation for the effective critical tension can then be written 



as 



= tiexp \{U(z max ) - U(z min )} /k B T] (7) 

where the potential is calculated with the tension r e . If the velocity is very small, the 
critical tension r e is very small and the escape time can be expanded in powers of the tension. 
For a vanishing tension, the escape time is t° e = tiexplUo/ksT) where t\ ~ (b/r (assuming 
that the potential has only one energy scale Uq ~ rb and one length scale b ~ Z\). The 
effective critical tension is then : 



kut° p 1 



N l + a(Tb/k B T)(t°Jt s (u)) 



(8) 



where t s (u) is the sticking time introduced in the previous section and a a numerical 
constant. This result is valid as long as the thermal escape time is small enough t° e <C 
^ft s (u) leading to £ « £ ~ (Na/b) 2 exp -(U /k B T). 

We assume for simplicity that the thermal escape time is larger than the Rouse time 
of the whole chain. The tension is thus uniform along the whole chain when the escape 
from the trap arises; the velocity u\ is then smaller than u c . In this limit of vanishing 
velocities, if the traps are independents, the friction force is Ff r = ^J- and vanishes as u 2 . 
If the escape time is smaller than the Rouse time only a fraction of the chain containing 
s ~ {kt Q e /C) 1 ^ 2 monomers is stretched. The critical tension is r e = {k^t^Y^u and the friction 
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force is Ff r — it also vanishes as u 2 . 

If the velocity is larger than Ui, the effective tension increases between a value of order 
ksT/b and r. The variation of the effective critical tension explicitly depends on the shape 
of the potential and cannot be expressed in simple analytical form. A simple expression can 
only be given if the velocity is large enough. In this case, the two positions z max and z m i n 
are close to the potential inflection point z\. One can then expand the potential around 

Z — Z\ 

U{z) =U 1 + (t- r e ){z - Zl ) - lyrzt i^)" (9) 

where 7 is a number of order one. The effective tension can be calculated from this 
expansion : 



^2/3 



Hi-iiog 2/3 (^)?(i-^) 1/2 A?)] (io) 



T T 



where we define the dimensionless number (3 = (4tzi) / (SksT^ 1 / 2 ) and the microscopic 
time t® = (£,2i)/(27 1 / 2 t). This expansion is valid if (3 is larger than log(t s /ti). Qualitatively 
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the effective tension is equal to the maximal force exerted by the potential if the sticking 
time t s is smaller than a time of the order of the thermal escape time. The effective tension 
increases smoothly (with a logarithmic dependence in u) between ksT/b and r. At first 
order we can use the power law found in the previous section in the framework of the zero 
temperature theory. 

A similar calculation can be made at larger velocities u 3> u c . If the thermal escape time 
is larger than the Rouse time, the effective critical tension is close to r and a law similar to 



eq. [10] is found (only the argument of the logarithm is slightly different). One can thus also 



use the zero temperature result for the friction force. 

Finally, if the traps are not independent one should compare the thermal escape time to 
the sticking time that, in this case, is of order d/u. The escape of the trap is dominated by 
thermal fluctuations if the thermal escape time is smaller. 



III. CONCLUDING REMARKS 

In this paper, we have studied a very specific model for the friction of a polymer chain on 
a solid surface. We have considered a two-dimensional gaussian polymer chain interacting 
with traps on a solid substrate. All the dissipation in this model is due to viscous friction 
in the solvent surrounding the chain and there is no dissipation associated to the substrate. 
Our main result is that, neglecting thermal fluctuations, the traps exert a finite friction 
force on the chain that does not vanish when the velocity goes to zero. This friction force 
originates in the elastic energy stored in the chain which is then dissipated by viscous friction 
in the liquid. At higher velocity the friction force decreases with the velocity due to the 
fact that the portion of chain that is stretched gets smaller. If thermal fluctuations are 
taken into account, the friction force formally vanishes at zero velocity. Nevertheless, for 
most practical cases the traps are strong enough and the effect of thermal fluctuations will 
appear only at infinitely small velocities. However, a very weak logarithmic dependence of 
the friction force on the velocity still remains, due to the fact that the trapping cross section 
depends logarithmically on the velocity. 

For a macroscopic sample consisting of a statistical ensemble of independent chains, the 
static friction coefficient is proportional to the average number of trapped chains. Further- 
more, our model gives an interesting result concerning the dependence of this coefficient 
with the sample history. During a friction experiment at velocity u, the average fraction of 
trapped chains is given by / ~ fs W M _ p or u <^ Uc this fraction is constant and equal to 
fo ~ tj • For u 3> u c it decreases with the velocity as f u ~ tJ— ■ In this case, if the motion 
is stopped fast enough, the static friction coefficient is lower than /q and increases with time 
until it reaches fo- 

Our model may look very simplistic, but we believe that many of the assumptions Ilmade 
here do not have a strong influence on the main results although they may change the 
quantitative value of the friction force. For example, we have assumed that only the last 
monomer of the chain (the chain end that is not pulled) interacts with the traps. It is 
easily checked that identical results are obtained (within a factor 2) if any monomer not 
too close to the chain end interact with the traps. We have explicitly performed the Rouse 
calculations in this case, they are slightly more complicated and we do not present them in 
this paper. Another extension that we would like to make in the future is the case where 
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several monomers of the same chain interact with the traps. There is of course a possibility 
that at the same time several monomers are trapped; this then leads to correlated release 
(cascades). We believe however that the same qualitative results would be obtained and that 
the friction force would remain finite at low velocities. We also have assumed a periodic array 
of traps. A more realistic approach would be to introduce disorder both in the positions of 
the traps and in the distance between traps. A weak disorder is clearly not important. If 
the traps are independent, then the averaging of the dissipation is straightforward. If the 
average distance between traps is small, the tension of the end monomer remains roughly 
constant as the tension fluctuations due to the trappings do not propagate to the pulled 
monomer. Disorder could play an important role if it introduces strong fluctuations (such 
as rare trapping events with a strong critical tension). It would then be interesting to study 
the fluctuations in the friction force. 

A last limitation is that we have only studied the stationary friction state. Transient 
effects are obviously very important and must be studied independently. 

It is also important to discuss the orders of magnitude of the predicted effects. Our 
more quantitative discussion of the release from the trap clearly shows that the solid-like 
friction can only be observed if the energy barrier is large rb/ksT ^> 1 otherwise the whole 
escape process is thermally activated and there is no solid friction. This could be achieved 
for example with a trap size of the order of a few monomer sizes b ~ 10a and rfc/Zc^T ~ 10. 
The critical velocity for a chain of 1000 monomers is then of the order of 0.1mm/ s i.e. in 
the accessible range for atomic force microscopy. One must however notice that the chain 
tensions can then be rather large and that it is not clear that the chain elasticity remains 
gaussian. This would of course change the quantitative value of the friction force but not 
the qualitative physics. 

In the future, we would like to use our very simple model to study the solid friction of a 
polymer gel moving on a solid substrate in order to give a quantitative interpretation of the 
recent experiments of Gong and Osada M 
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APPENDIX A: TRAPPED CHAIN. 



It is convenient to perform all calculations using the Laplace transformation, denned by : 



f + OO 

z(s,p) = / z(s,t)e~ pt dt 
Jo 



d 2 



For the trapped chain, the following Rouse equation has to be solved: 

z free 

with the boundary conditions : 

z(s = 0,p) = 

z(s = N,p) = — 
p A 



The solution of this equation is given by: 

sinh \/T R p 



Zd(s,p) 



V 



1 + \pT» ( 



+ 



(s-N) 
N 



z d (s,t) = uT R 



t S 1 / a \ - t—^ 

T R ~N + 2\Nl ^ 



sinh ^jT^p 
- 2 +2? (-1)' 

7T° 



n=l 



TV 



1-e 



-n 2 7rV 

1 i 



n 



' R ) sin [ n ^( s - N ) 



and the tension along the chain is given by: 



T c h(s,t) = k 



uT, 



R 



N 



t S 2 ±2? 1 / -n 2 7r 2 -*-\ ( n , 



At large time t ^> Tr, the tension on the trapped monomer behaves as : 

T tr {s = 0,t) ~ -^Ut 

which corresponds to the tension of N springs under the extension ut. 
For short times t -C Tr, the tension on the trapped monomer behaves as 



T tr (s = 0,t) 



2 k 



uT R \lJ- 



T, 



R 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



As given in the main text, the sticking time is obtained when this tension is equal to the 
critical value r. 



In this case, the number of extended springs is s (t) ~ Ny 



tr- 
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APPENDIX B: DETACHMENT 



After the detachment of the trapped monomer, the chain relaxes to the steady state, and 
its elastic energy is dissipated. The Rouse equation becomes: 



where z d (s) ~ uT R 



d 2 

C (pz{s,p) - z d (s)) = k—z(s,p) 



(Bl) 



(i n erfc(x) is a repeated inte- 



gral of the error function defined in pB|) gives the conformation of the chain at the time 
where the chain detaches form the trap . The boundary conditions are: 



dz 
ds 



8 = 0,p) = 



Z[S 



z{s = N,p)=uT R (- + -?-)- + - 
\2 T R J p p 2 

The solution of ( pi|) is given by: 

i It / r s 1 

,p) = Au + {s,p) + Bu-(s,p) — J— (m_(s,p) / z^s'y^*^ ds' - u + {s,p) 

2JM V V V jo 



(B2) 



z d (s')e~^^ds' 
(B3) 



where u± = e ± ^ pTR n , A and B are two constants determined by ( |B2|) . The first boundary 
condition gives A = B = ^. In the limit of large times, the second condition gives a ~ 



mi 



cosh(v / p7jj) 



. At large times, the tension on the pulled monomer is given by: 



u 



f d (N,p) ~ AC- + 8i s erfc(0)N(uJ-?-t 



V 



T, 



R 



r d {N, t) ~ A(n + 8i 3 erfc(0)N(uJ ^-t a S(t) 



T, 



(B4) 



By integrating this tension over time one finds the dissipated energy given in the main text 
(equation |j) . 
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APPENDIX C: HYDRODYNAMIC INTERACTIONS. 



The calculations presented in the main text of this paper have been made using the 
Rouse dynamics for the polymer chains on the solid surface. This ignores two effects : the 
direct friction of the chains on the solid surface and the hydrodynamic interactions between 
monomers. The hydrodynamic interactions close to a surface are difficult to calculate explic- 
itly. However if we remain at a scaling level, it is sufficient to introduce the local stretching 
blobs [^7[. If the local tension of the chain is r c h, the local blob has a size £ t ~ k B T /r c h and 



contains g t ~ (£ t /a) 2 monomers. The friction force on a blob is Qin]C, t and the hydrodynamic 
interactions are screened at the scale of the blob so that the dynamics of the chain of blobs 
can by studied using the Rouse model; the total friction is the sum of the frictions over all 
the blobs. As in the Rouse model, a free chain pulled by one end with a velocity u has on 
average a trumpet like shape studied in details by Brochard and her coworkers |2~8 |. The 
Zimm equation of motion is obtained by writing a force balance on the blobs 

z dz 

k d? = 6nr]U dS (C1) 

The average position of monomer s is then 

fnrniiQ 

(C2) 



Zfree \S ) 



Qnrju 



1/2 f 6nr]us 
I exp( - ) - 1 



The force balance on a section of chain containing s monomers close to the free end point 
is obtained by integration of the equation of motion = 6tt7]uz. Note that the tension 
increases very rapidly from the free end point and that the nonlinear elasticity certainly 
becomes relevant. We do not discuss this effect here and refer to reference |28| . 



For a chain trapped by one end on the surface, the sticking time does not depend on 
the dynamical model at low velocities (if it is larger than the characteristic relaxation time 
of the chain which is the Rouse time of the chain of blobs) t s = Nr/(uk). The dissipated 
energy is still the elastic energy of a chain with uniform tension r. 

At high velocities, the sticking time is obtained by a scaling argument similar to that 
used in the Rouse model. When the end point escapes from the trap after a time t s , a chain 
section of s(t s ) monomers is stretched with a tension r. The characteristic relaxation time 
of this chain section is the Rouse time of the corresponding chain of blobs t s = 67r ^ T ■ Using 
the integrated equation of motion, we find the sticking time t s = §J 2 ■ It nas the same 
dependence on the velocity as in the Rouse model but a different dependence on the tension 
which is due to the fact that one must now consider the chain as a Rouse chain of Pincus 
blobs and not of monomers. 

The crossover velocity is obtained by comparing the low and high velocity results u c = 
6 J° N and the energy dissipated on a trap at high velocity is the elastic energy of a chain 

section containing s(t s ) monomers: Wd = Here also we find the same dependence on 

velocity but a different tension dependence if we compare to the Rouse model. 
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APPENDIX D: TRAPPING CROSS SECTION 



In this appendix, we study the trapping cross section for the free chain end by the 
potential well of a trap on the surface. We focus on the undisturbed part of chain end of 
size dr. At very low velocity u <C Mj = R/Tr (Tr is the Rouse time) and 5r is the chain 
radius R. At higher velocity the chain is stretched by the hydrodynamic friction and Sr is 
equal to the local stretching blob size (Sr ~ w -1 / 3 ). If we take the origin at the center of the 
first blob containing the free chain end, the surface monomer concentration at a distance r 
is : 

c(r) ~ Sr~ 2 exp — (r/Sr) 2 (Dl) 

The trap acts as an absorbing site and we consider that the chain end gets trapped when 
it diffuses to the surface of the trap of radius b (figure [l]). If the trap is at a distance r from 
the average position of the end point, the concentration gradient close to the trap is c(r)/b. 
The total flux of end point towards the trap is thus given by : 

q = 2irbDc(r)/b (D2) 

The relevant diffusion constant D is equal to the Rouse diffusion constant of the whole 
chain D = ksT/^NQ at very low velocity and to the Rouse diffusion coefficient of the 
stretching blob D = k^T j (sQ (s is the number of monomers in the blob) at high velocity. 
Note that in two dimensions there are logarithmic corrections to this result as there is 
no stationary diffusion state. If the impact parameter (the minimum distance between 
the trap and the average chain end position) is c and if the distance is minimum at time 
t = the distance r varies with time as r 2 = c 2 + u 2 t 2 . The probability p(c) that the 
chain with an impact parameter c gets trapped is obtained by integration over time of the 
flux q, p(c) ~ m exp — (c/R) 2 . This probability is of course smaller than one. The trapping 

probability is thus equal to one if c < R log 5 (uq/u) and decays as a gaussian for larger values 
of c. The trapping cross section (a length in two dimensions) is obtained by integration over 
the impact parameter 

a ~ R\og^(u Q /u) (D3) 

At higher velocity u > uq, the trapping probability is always smaller than 1 and one 
can integrate the flux q to find the trapping probability and then the cross-section. One 
must however take into account that 5r is now the size of the stretching blob and that the 
diffusion constant is not the diffusion constant of the whole chain but the one of a chain of 
size 5r. The result is 

The average distance between trapping events can then be calculated from the surface 
density of traps n, d = 1/an. 
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Figures caption 



Figure 1: Polymer chain pulled at constant velocity u on a periodic 2D lattice of 
traps. 

Figure 2: Force (normalised by N^w) applied on the trapped monomer (s=0, black) 
and on the pulled one (s=N, grey) versus time (normalised by Tr), according to 
equation [| 

Figure 3: Friction force in a (u, d) diagram for the many traps problem at zero 
temperature. 

Figure 4: Trapping potential a) V(z) and b) U(z) = V(z) — r tr z. 
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FIGURES 
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FIG. 2. 
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FIG. 3. 




